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ABSTRACT 

We present NuSTAR observations of the powerful radio galaxy Cygnus A, focusing on the central 
absorbed active galactic nucleus (AGN). Cygnus A is embedded in a cool-core galaxy cluster, and 
hence we also examine archival XMM-Newton data to facilitate the decomposition of the spectrum into 
the AGN and intracluster medium (ICM) components. NuSTAR. gives a source-dominated spectrum 
of the AGN out to > 70keV. In gross terms, the NuSTAR spectrum of the AGN has the form 
of a power-law (F ~ 1.6 — 1.7) absorbed by a neutral column density of Njj ~ 1.6 x 10 23 cm -2 . 
However, we also detect curvature in the hard (> lOkeV) spectrum resulting from reflection by 
Compton-thick matter out of our line-of-siglit to the X-ray source. Compton reflection, possibly from 
the outer accretion disk or obscuring torus, is required even permitting a high-energy cutoff in the 
continuum source; the limit on the cutoff energy is E cut > lllkeV (90% confidence). Interestingly, 
the absorbed power-law plus reflection model leaves residuals suggesting the absorption/emission 
from a fast (15,000 — 26,000kms -1 ), high column-density {Nw > 3 x 10 23 cm -2 ), highly ionized 
(£ ~ 2, 500 erg cm s^ 1 ) wind. A second, even faster ionized wind component is also suggested by 
these data. We show that the ionized wind likely carries a significant mass and momentum flux, 
and may carry sufficient kinetic energy to exercise feedback on the host galaxy. If confirmed, the 
simultaneous presence of a strong wind and powerful jets in Cygnus A demonstrates that feedback 
from radio-jets and sub-relativistic winds are not mutually exclusive phases of AGN activity but can 
occur simultaneously. 

Subject headings: accretion, accretion disks — galaxies: clusters: intracluster medium — galaxies: 
jets — X-rays: individual: Cygnus A 


1. INTRODUCTION 

Ever since its discovery in the early days of ra- 
dio astronomy, the po werful radio galaxy Cygnus A 
(lHargrave fc Rvld I1974H has been an important proving 
ground for our models of active galactic nuclei (AGN) 
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and the relativistic jets that they produce. It is extremely 
radio luminous, with a 178 MHz luminosity more than an 
order of magnitude larger th an any other 3C sourc e in 
the local (z < 0.1) universe (iCarilli fe Barthellll996D - 
indeed, one needs to reach out to the z ~ 1 universe 
before finding many other sources of comparable radio 
luminosity. Much of the radio emission originates from 
two edge-brightened radio lobes that are fed by relativis¬ 
tic and highly collimated back-to-back jets produced by 
the AGN, making Cygnus A th e arc hetypal exam ple of 
a Fanaroff-Riley Type-II (FRII; 'Fanaroff & R.ilevill97 4) 
classical double radio galaxy. 

Studies of the central engine of Cygnus A are hampered 
by heavy obscuration. The di scovery of a broad MgH line 
in the ultraviolet spectrum (lAntonucci. Hurt fc Kinnevl 
I1994T) a nd broad Ha in the polarized optical spectrum 
(lOgle et al.l Il997t) makes it clear that the central engine 
of Cygnu s A is a broad-line qua sar, as previously sug¬ 
gested bv lDiorgovski et al . (199 11 on the basis of infrared 
(IR) imaging. Interestingly, Cygnus A appears to be only 
a modestly powerful q uasar despite powering a n ultra- 
luminous radio-source dBarthel fc Arnaud 199611 . A re¬ 
cent demonstration of this was given bv iPrivon et al.1 
( 2012 ) who model the radio-to-IR spectral energy dis¬ 
tribution (SED), concluding that the IR luminosity is 
dominated by the AGN with a bolometric luminosity of 
Lboi « 4 x 10 45 ergs _1 , typical of a low-power quasar. 
For a black hole mass of M = (2.5 ± 0.7) x 10 9 M Q , 
determined via Hubbl e Space Telescope spe ctroscopy of 
a nuclear gas disk (|Tadhunter et al.l 12003( 1 . this bolo¬ 
metric luminosity implies a modest Eddington ratio of 
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L/L^dd ~ 0.01. Substantially more power is believed to 
emerge in a kinetic form associated with the rela t ivistic 
jets; using the dynamics of the cocoon, llto et alJ (|2008H 
estimate a total jet power in the range Lj ~ (0.7 — 4) x 
10 46 ergs _1 . An additional reason for the anomalous 
radio-loudness of Cygnus A is its unusual environment 
it is host ed by the cD galaxy of a cooling core cluster of 
galaxies (lArnaud et all 1984 iRevnolds fe Fabianl [l996t 


iBarthel fc Arnaudl 11990 : iSmith et al.l l2002fl~ The high 

pressure core of the intracluster medium (ICM) provides 
a working surface for the jets and confines the radio lobes, 
thereby greatly increasing the synchrotron emissivity of 
the shocked jet plasma. 

Hard X-ray spectroscopy is a powerful way to probe the 
central engines of deeply buried AGN such as Cygnus A, 
although the presence of the X-ray luminous ICM com¬ 
plicates such studies. While the X-ray emission from 
the Cygnu s A cluster was dete cted by the Einstein Ob¬ 
servatory (lArnaud et alJ ITOSlI) ■ the harder spectral re¬ 
sponse of the medium-energy (ME) telescope on EX- 
OS AT was required to first detect the additional hard 
X-ray emission from the absorbed nucleus (jArnaud et all 
1987). Ginga permitted a robust detection of the AGN 
continuum out to 20keV (llleno et al.l 119911k the contin¬ 
uum was well described by a power-law with photon in¬ 
dex r = 2.0 ± 0.2 and 2-10 keV intrinsic luminosity of 
1 x 10 45 ergs _1 absorbed by a neutral column density 
of Nh = (3.8 ± 0.7) x 10 23 cm -2 . This basic picture, 
albeit with a somew hat flatter photon index (T « 1.5), 
was confirmed bv I Young et all ( 2002 .) who used a combi¬ 
nation of Chandra and the Rossi X-ray Timing Explorer 
(RXTE) to map the AGN emission out to 100 keV. By 
using its superior spatial resolution to isolate the nuclear 
emission, Chandra found X-ray reprocessing of the AGN 
continuum in the form of a weak 6.4 keV fluorescent line 
of cold iron. However, the non-imaging nature of the in¬ 
struments on RXTE meant that the hard X-ray spectra 
were strongly background dominated. 

In this paper, we present observations of Cygnus A by 
the Nuclear Spectros copic Telescopic Array ( NuSTAR ; 
lHarrison et al.l 12013 :1. The imaging capability of this 
focusing X-ray obser vatory , with a half-power diame¬ 
ter of just 60arcsec dMadsen et al.l [20T5h . allows us to 
produce a high signal-to-noise (S/N), source-dominated 
spectrum of Cygnus A out to almost 80keV. The un¬ 
precedented quality of this spectrum allows us to search 
for the presence of Compton reflection signatures, con¬ 
strain any high-energy cutoff of the continuum, and 
search for highly ionized outflows for the first time in 
this keystone object. 

An important complication, however, is the contribu¬ 
tion of the ICM emission in the softest bands of NuS¬ 
TAR and the all important iron K-shell band. Thus, 
our NuSTAR analysis must be informed by additional 
soft X-ray (0.5-10 keV) imaging-spectroscopy that allows 
us to construct a spectral model for the ICM. At first 
glance, the moderately deep Chandra observ ations (to¬ 
taling 200 ks) of Cygnus-A (Smith et al.ll2002 1 would ap¬ 
pear to be the datasets of choice for this exercise. How¬ 
ever, the AGN itself suffers extreme photon pile-up dur¬ 
ing these observations (compromising the integrity of any 
global ICM+AGN spectrum), and the hard X-ray wings 
of the Chandra PSF noticeably contaminate the inner 



Figure 1 . Folded background-subtracted spectra from the FPMA 
(red filled circles) and FPMB (blue filled circles), along with the 
background spectra (FPMA background shown as red open circles 
and FPMB as blue open circles). 


10 arcsec of the ICM with AGN emission (compromising 
the ability to spatially isolate the ICM emission). Thus, 
in order to obtain the best constraint on the ICM emis¬ 
sion contaminating our NuSTAR view of the nucleus, we 
turn to archival XMM-Newton observations. 

This paper is organized as follows. The NuSTAR and 
XMM-Newton observations and the basic data reduction 
steps are described in Section 2. Section 3 gives a pre¬ 
liminary discussion of the NuSTAR image of Cygnus A, 
although we defer a detailed imaging analysis, includ¬ 
ing the search for extended hard X-ray emission from 
non-thermal particles, to a future publication. Section 4 
presents our analysis of the XMM-Newton data. With 
the XMM-Newton ICM model in hand, we describe the 
analysis of the NuSTAR spectrum in Section 5. Section 6 
discusses the astrophysical implications of our results. 
Throughout this work, we ass ume the stand ard Planck 
cosmology (iPlanck Collaboration et ahll2013f ). At a red- 
shift of z = 0.056, this places Cygnus A at a luminosity 
distance of 237 Mpc, with a linear-angular conversion of 
1.0 kpc per arcsec. 

2. OBSERVATIONS AND DATA REDUCTION 
2.1. NuSTAR 

NuSTAR observed Cygnus A on 17/18-Feb-2013 and 
l-Mar-2013 with a total observation length of 80 ks and 
40 ks respectively. The data from both focal plane mod¬ 
ules (FPMA and FPMB) and both observations were 
reprocessed and cleaned using the most recent version 
of the NuSTAR pipeline (within HEASOFTvl6.6) and 
calibration files (CALDB version 20140715), resulting in 
43.6ks and 20.7ks of good on-source exposure for the 
first and second observation, respectively. 

Although we briefly discuss the imaging data in the 
next section, this paper focuses on the AGN at the heart 
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of Cygnus A. We extracted spectra and light curves from 
a 30 arcsec radius circular region centered on the nucleus 
(dashed circle in Fig. [2]). Our chosen extraction region 
is somewhat smaller than usual in order to minimize the 
contribution of the ICM. Background spectra and light 
curves were extracted from two circular regions that are 
free of any obvious point sources and flank the nucleus by 
4arcmin. To remain within the well-calibrated regime, 
this paper considers the 3-70 keV NuSTAR spectrum and 
we rebin the data to a minimum of 20 photons per en¬ 
ergy channel in order to facilitate the use of x 2 -statistics. 
The (folded) background-subtracted source spectra for 
the two FPM are shown in Fig. [T] along with the corre¬ 
sponding background spectra. We see that the spectra 
are source-dominated across the whole 3-70 keV — the 
background contributes only 5% at 7keV increasing to 
25% at 70keV. 

Examination of the background-subtracted light curves 
reveals no evidence for time variability within an obser¬ 
vation; this has been tested with lightcurves employ¬ 
ing 500 s, 1000 s, and 2000 s bins. This result stands, 
even when we restrict our attention to energies above 
10 keV (thereby largely removing the contribution of the 
constant ICM component). Additionally, between the 
two observations, the average background-subtracted 3- 
70 keV count rate of each FPM is found to be constant 
within errors. Hence, we combine the data from the 
two observations to produce a single spectrum (for each 
FPM) with a total exposure time of 64.3 ks. There is 
a 1.7% offset in the count rate between the two FPMs 
(with count rates of 0.668 ± 0.003 cps for FPMA and 
0.656 ± 0.003 cps for the FPMB) that we attribute to 
residual errors in the flux calibration of the two instru¬ 
ments. In all of the spectral fitting performed in this pa¬ 
per, we fit for the cross-normalization of the two FPMs 
thereby removing this cross-calibration error. 

2.2. XMM-Newton 

XMM-Newton observed Cygnus A on 14-0ct-2005 
(22.5ks) and 16-0ct-2005 (18.8ks). Here we consider 
just the imaging spectroscopy data from the EPIC-pn de¬ 
tector. These data were extracted from the HEASARC 
archives and reprocessed using the epchain tool within 
SASv 13.5.0. After applying standard filtering criteria 
(described in the XMM-Newton ABC guide), the first 
and second observations yielded 18.3 ks and 15.0 ks of 
good on-source data, respectively. 

Our motivation for examining the EPIC-pn data is to 
give context to our NuSTAR study. Hence, for each of 
the two observations, we extract spectra from the NuS¬ 
TAR extraction radius, a 30 arcsec radius circular region 
centered on the nucleus. Background spectra were ex¬ 
tracted from a region approximately 4 arcmin from the 
nucleus that is free from any obvious point sources. The 
source and background spectra were appropriately scaled 
(using the backscale tool), and response matrices and 
effective area curves generated (using the rmfgen and 
arfgen tools). Finally, given the proximity in time of 
the two observations and the lack of any detected vari¬ 
ability between these two observations, the two spectra 
were combined into a single spectrum with an exposure of 
33.3 ks. To remain within the well-calibrated regime, this 
paper considers the 0.5-10 keV EPIC-pn spectrum and 
we rebin the data to a minimum of 20 photons per en¬ 


ergy channel in order to facilitate the use of x 2 -statistics. 

3. IMAGES 

A detailed study of the hard X-ray imaging data for 
Cygnus A, including constraints on non-tliermal X-ray 
emission from the Cygnus A cluster, is deferred to a fu¬ 
ture paper. Here, we present a preliminary discussion 
of the imaging data to the extent required to set up the 
analysis of the nuclear spectrum. Figure [2] shows the 0.5- 
lOkeV XMM-Newton/ EPIC-pn image of the Cygnus A 
along with NuSTAR images made in four bands (3-5 keV, 
5-10keV, 10-20keV and 20-40keV). For reference, we 
also show the 0.5-5 keV Chandra /ACIS imag43- All 
frames are overlaid with contou rs of the classic 6cm Very 
Large Array (VLA) map of iPerlev. Dreher fe Cowanl 
(11984H as well as the NuSTAR extraction radius (dashed 
circle). Each frame has been scaled with a logarithmic 
color map spanning a dynamic range of 100. 

The Chandra image clearly shows the prolate el¬ 
lipsoidal cocoon b lown by the radio galaxy activity 
dWilson. Smith fc Youndi2006f ). The emission visible be¬ 
yond this cocoon is the ambient ICM. As previously 
noted by iSmith et al.l (12002H , the cocoon appears to be 
wrapped with X-ray bright filaments, the origin of which 
is unclear. We see that the NuSTAR extraction ra¬ 
dius fits within the radio lobes and encompasses the 
inner portions of the jet-blown cocoon. In the XMM- 
iVeiafon/EPIC-pn image, the cocoon structure is less dis¬ 
tinct but the X-ray emission from the radio hot-spots is 
easily visible. In NuSTAR, as expected, the extended 
cluster emission is also prominent in the 3-5 keV image. 
However, as can be see from the bottom panels of Fig. [2J 
the nuclear point source becomes increasingly dominant 
as one considers harder and harder bands 1 ”!. As we will 
see, this is completely in line with our spectral analysis. 

4. XMM-NEWTON SPECTRAL CONSTRAINTS ON THE 
AGN AND INNER ICM 

We now present an analysis of the X'MM- 
Newton/EPIC-pn spectrum of Cygnus A. A spatially 
resolved investigation of the interaction of the radio 
lobes with th e ICM using these data h as been presented 
previously bv iBelsole fc Fabianl (12007I 1. However, to the 
best of our knowledge, the study presented here is the 
first published analysis of the XMM-Newton view of 
the Cygnus A nucle u s itself. Based on pr e vious works 
(lYoung et al.l 120021 lArnaud et al.l 119871 : lUeno et all 
l1994h . we expect the X-ray spectrum within our extrac¬ 
tion radius to be the superposition of thermal emission 
from the optically thin plasma of the ICM and a heavily 
absorbed AGN spectrum (all abso rbed by a Galactic 
column of Nh = 3.5 x 10 21 cm -2 ; iDickev &; Lockmanl 
fl990h . As we will see, the EPIC-pn derived constraints 
on the parameters of the AGN are rather poor due 
to the combination of the heavy absorption and the 
truncation of the EPIC bandpass at lOkeV. However, 
our main motivation here is to construct an ICM model 

15 This image is formed by merging the level-2 events files from 
ObsIDs 5830, 5831, 6225, 6226, 6228, 6229, 6250, and 6252. 

16 We note that there is only a w eak energ y depe nden ce to the 
point spread function of NuSTAR IMadsen et all [201511 . with a 
half-power diameter increasing by approximately 10% at the softest 
energies (< 4.5 keV). 
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Figure 2. X-ray images of the Cygnus A field overlaid with contours of the lPerlev. Dreher fe Cowanl (i 19841 ') 6 cm radio map (cyan) and 
our XMM-Newton/NuSTAR spectral extraction region (dashed circle of 30 ,/ radius). Each frame has been scaled with a logarithmic color 
map spanning a dynamic range of 100. 


that can be taken over into our NuSTAR analysis of the 
AGN. 

We find that the IC M emission is w ell described by 
two apec components (ISmith et al.ll200ll ) with indepen¬ 
dent temperatures (Tj, T 2 ). Driven by physical consider¬ 
ations, we assume that the two ICM components have a 
common metal abundance (Z). In this paper, all abun¬ 
dances are re ferenc ed to the sola r abu nd anc e (Z Q ) set 
described by iWilms. Allen fc McCravl (120001) . For all 
spectral models described below, the inclusion of a third 
temperature ICM component fails to give a significant 
improvement in the goodness of fit. 

We begin by approximating the AGN emission with 
a simple power-law (PL) continuum absorbed by a neu¬ 
tral/cold absorber at the redshift of Cygnus A (cABS) 
with column density Nh and solar abundances. Both 
the intrinsic absorber and the Galactic absorbing column 
(which also affects the ICM emissio n) are descr ibing us¬ 
ing the tbabs code of IWilms. Allen fe McCravl (|2000f ). 
This model provides a good description of the broad¬ 
band spectrum (model ICM+cABS[PL] in Tablejl]). The 
residuals do indicate, however, a strong unmodeled emis¬ 
sion line between 6-7keV (rest-frame). Adding a Gaus¬ 
sian emission line (emLINE) into the spectral model 
leads to a dramatic improvement in the goodness of 
fit (Ay 2 = —243 for three additional model parame¬ 
ters). The line is narrow (er < 70 eV corresponding 


to a FWHM< 7700kms -1 ), relatively strong (equiva¬ 
lent width W = 98 ± 12 eV), and has a centroid en¬ 
ergy precisely that expected from the fluorescent Ka 
line of cold iron (E = 6.39 ± 0.01 keV). This model 
(ICM+cABS[PL+emLINE]) leaves no gross residuals in 
the fit. 

In their Chandra study of Cygnus A, lYoung et al.1 
(12002D find soft power-law emission from a bipolar source 
spanning the nucleus that they identity as scattered AGN 
emission. Hence, we must explore whether the ICM pa¬ 
rameters can be skewed by the addition of such a scat¬ 
tered component. Adding a power-law component that 
is not obscured by the AGN absorber with a photon in¬ 
dex that is tied to that of the main AGN continuum 
b ut with a normalization decreased by a factor of 100 
(lYoung et all 1200^ 1 has essentially no impact on the fit 
or the ICM parameters. Allowing the normalization to 
be free results in a small and insignificant improvement 
in the goodness of fit (Ay 2 = —7 for 2 additional de¬ 
grees of freedom) and a best fitting ICM abundance of 
Zp e = 2.5^QgZQ. However, the normalization of this 
scattered component is then 25% of the primary con¬ 
tinuum — this is 25 x higher than found by Chandra , 
strongly suggesting that we should reject this possibility. 
We conclude that soft/scattered emission from the AGN 
does not strongly bias our ICM model. 

Physically, iron lines such as seen in our X'MM- 
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Observed energy (keV) 



Figure 3. Fit of the ICM+cABS[PL+REFL] model to the XMM-NewtonfEPlC-pn data. Left panel : 0.5-10 keV folded EPIC-pn spectrum 
decomposed into contributions from the cool (red) and hot (brown) ICM components, plus the absorbed and reflected AGN power-law 
component (blue). Right panel : Zoom-in on the iron-K region of the folded spectrum. The resolution of the EPIC-pn allows contributions 
to the iron complex from the fluorescent line of cold iron (6.40 keV rest-frame; 6.06 keV observed) as well as the radiative recombination 
lines of iron-25 and iron-26 (6.67 keV and 6.97 keV rest-frame; 6.32keV and 6.60 keV observed) to be discerned. The cold iron line is modeled 
as AGN reflection, whereas both ionized iron lines can be attributed to the hot component of the ICM. 


Spectral Model 
ICM-bcABS(PL) 


ICM+cABS(PL+emLINE) 


ICM+cABS(PL+REFL) 


Parameters 

kTi = 1.69+o!io keV> kT 2 = 6.74^“keV, Z ICM = 1.33+“;^Z 0 
N H = 41.2±|;® x 10 22 cm" 2 , T = 
kTl = 1.70^ keV, kT 2 = 7.04±“;« keV, Zjcm =''l.40t£}£Z Q 

AT h = 34.1±|;g x 10 22 cm“ 2 , T = 1.43 ± 0.11 
E liue = 6.39 ± 0.01 keV, o- Une < 70 eV, Wi ine = 98 ± 12 eV 
kTi — l-70^Q Qg keV, kT 2 = 6.70^Q - 3g keV, Zicm — 1.30 dz O.IO^q 
Nn = 31.2 +11 x 10 22 cm" 2 , T = 1.51±g 

n = igtojf. lo §^ < 2 - 15 ' Z L = z & 


X 2 /dof 

1712/1605 


1469/1602 


1461/1603 


Table 1 

Spectral fits to our XMM-Newton EPIC-pn spectrum. See text (0 for a detailed discussion of the spectral models. All energies are 
quoted in the rest-frame of Cygnus A (which has a cosmological redshift of z = 0.056). The ionization parameter is i n units of 
ergcms -1 . Superscript / denotes a fixed parameter. All errors are quoted at the 90% level for one interesting parameter. 


Newton spectrum are expected to arise fr om X-ray re¬ 
flection by Compton-t hick cold matter (|Baskol [L978S; 

I George fe Fabiaiil [1991); this matter can be identi¬ 
fied with the outer regions of the accretion disk or a 
Compton-thick core to the obscuring torus. Hence, we 
replace the Gaussian line component with a cold X-ray 
ref lection mo del as calculated by the xillver code of 
iGarcfa et all 1 2013 ) — we set the abundance to solar 
but allow the ionization state of the reflector and the 
reflection fraction to be free parameters. This model 
(ICM+cABS[PL-|-REFL]) provides a very comparable fit 
to the spectrum (Table [lj. The AGN parameters, espe¬ 
cially those describing the X-ray reflection, are poorly 
constrained. Allowing the metallicity of the reflector to 
be a free parameter fails to improve the fit and permits 
a strong degeneracy between the metallicity and the re¬ 
flection fraction. This is readily understood - in the 
XMM-Newton band, the principle diagnostic of reflec¬ 
tion is the energy and strength of the iron line, and a 
given line strength can be achieved by weak reflection 


from gas with high metallicity, or strong reflection from 
gas with low metallicity. 

Figure [3] (left) shows the folded EPIC-pn spectral fits 
decomposed into its model components together with the 
fit ratios. The ICM emission (consisting of kT = 1.7 keV 
and 6.7 keV components) dominates the spectrum below 
4keV, with the AGN emission becoming dominant at 
higher energies. The iron line complex (Figure |3] right) 
has distinct peaks corresponding to the 6.4keV fluores¬ 
cent line of neutral iron (from X-ray reflection within the 
AGN), and the radiative-recombination lines of helium- 
and hydrogen-like iron at 6.67 keV and 6.97keV respec¬ 
tively (from the hot ICM). 

5. ANALYSIS OF THE NUSTAR SPECTRUM 

5.1. Basic Characterization : The Detection of 
Compton Reflection 

For an initial impression of the NuSTAR spectrum, we 
compare the data with a spectral model consisting of a 
two-component ICM and a power law continuum mod- 
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Spectral Model 

Parameters 

X 2 /dof 

ICM+cABS(PL) 

N h = 25.8i2 2 x 1022 cm 2 > r = 1.57 ± 0.02 

1410/1148 

ICM+cABS(PL+emLINE) 

N h = 19.8+ 2 g x 10 22 cm” 2 ,T = 1.54 ± 0.02 

Eemis = 6.34l°;° 3 keV, o- emia = 0.01 keV*, W emia = 147±|| eV 

1232/1146 

ICM+cABS(PL+REFL) 

N h = 16.9+) ( x 10 22 cm- 2 , T = 1.77 ± 0.03 

log Crefi = 2.7l°;),77 = 0.50l°;J°,Z refi = 2.3l°; 3 

1237/1145 

ICM+cABS(PL+REFL+absLINE) 

N h = 22.4 ± 0.7 x 10 22 cm" 2 , T = 1.76l°;° 3 

logged < 1.1, n = 0.43 + 0.04, Z refl = 2.5l°;° 

1216/1143 


£ abs = 7.18 + 0.08 keV,<r{ bs = 0.01 keV, VF aba = -501)® eV 


ICM+cABS(PL+REFL+emLINE) 

N h = 19.6l) | x 10 22 cm" 2 ,T = 1.671q’q5 
logged = 2.7l°;i,7e = 0.4o!°;°£, Z refl = 0.5l° p 16 

1188/1143 


Semis = 6.30 ± 0.04 keV, <rf mi8 = 0.01 keV, IF emia = llolU eV 

Tvr - + 1.Y . . ... 9 -,^+U.U3 



ICM+cABS(PL+REFL+emLINE+absLINE) N H = 19.6+2 I x 1022 = l^l^g 1172/1141 


logged = 2.7±°- 1 3 ,nm 0.391°;°!, Z rea = 0.5l° p 21 
Semis = 6.30 ± 0.04 keV, af mis = 0.01 keV, lF emia = 107± 2 ^ eV 
E aba = 7.24 ± 0.10keV,<rf bs = 0.01 keV, W aba = -35±^ eV 

ICM+cABS*wABS(PL+wEMIS+REFL) N H = 17.6 ± 2.3 x 10 22 cm" 2 ,F = 1.60lgg) 1163/1140 

Nwabs > 3.4 x 10 23 cm -2 , log= 3.4l°;®, D W abs = -1.91J ;° x 1q4 kms_1 
Nwemis = 2.9l°;| x 10 23 cm -2 , wwemis = +3.2 ± 0.3 x 10 4 kms _1 

_ loggrefl = 2.7l°; 3 ,7e = 0.36l°;)®,Z refi <13 _ 

ICM+cABS*2wABS(PL+wEMIS+REFL) N H = 16.9 ± 1.4 x 10 22 cm" 2 , f = 1.47lgog 1144/1137 

N W absi > 3.0 x 10 23 cm -2 , log £vvi = 3.4l°J, u Waba i = -2.0l°jj x 10 4 kms _1 

Nwabs 2 > 3.6 x 10 22 cm -2 , log S^w 2 = 3.01J ;% abs 2 = -5.0l 3 ;° x 10 4 kms _1 

Nwemis = 5.0l)®g° x 10 23 cm -2 , t>wemis = +3.1 ± 0.2 X 10 4 kms -1 
loggrefl = sstl ln = onlgjj, z refl < 3.0 


Table 2 

Spectral fits to the 3-70 keV NuSTAR data. See text (EH and E3 for a detailed discussion of the spectral models. All energies are 
quoted in the rest-frame of Cygnus A (which has a cosmological redshift of z = 0.056). All ionization parameters (£+/ • £yyi • /u/ 2 , £ r ed) are 
in units of ergcms -1 . Superscript / denotes a fixed parameter, and sub-script p denotes a parameter that has hit the upper/low limit of 
its allowed range. All errors are quoted at the 90% level for one interesting parameter. 



Figure 4. Ratio of the NuSTAR FPMA (red) and FRMB (blue) 
data to a spectral model in which the AGN is described by a simple 
absorbed power-law (model ICM+cABS[PL] in Tablc[2]). The spec¬ 
tral model includes a two-temperature ICM component with shape 
(but not normalization) fixed to that found in the XMM-Newton 
analysis. 

ified by intrinsic cold absorption. The “shape parame¬ 
ters” of the ICM spectrum (i.e. the temperatures, the 
common abundance, and the relative normalizations of 


the two components) are fixed to those of our best-fitting 
model for the XMM-Newton data. The total normaliza¬ 
tion of the ICM emission is allowed to depart from the 
XMM-Newton value to account for the different point 
spread functions and any possible flux cross-calibration 
issues. This model is not a particularly good description 
of the data [model ICM+cABS(PL) in Tabled- As can 
be seen in Fig. Q] (also Fig. [5^), this model leaves obvi¬ 
ous residuals in the form of an emission line at iron-Ka 
energies and a weak broad hump peaking at 20-30 keV. 
This high-energy hump is suggestive of Compton reflec¬ 
tion from optically thick matter. 

Adding a narrow Gaussian emission line to the model 
leads to a dramatic improvement in the goodness of 
fit (Ay; 2 = —178 for two additional model parameters; 
model ICM+cABS[PL+emLINE] in Table [2]). As ex¬ 
pected, the putative Compton reflection hump is still not 
captured by this model, leading to the subtle high-energy 
hump visible in panel-b of Fig. [5] Of more interest is the 
fact that the simple emission line model leaves significant 
unmodeled complexities in the iron-K band (6-8keV), 
and that the best-fitting energy of the line (rest-frame 
E = 6.34to(o4 keV at the 90% confidence level) is in¬ 
consistent with that expected from any charge-state of 
iron. While it may be tempting to interpret this as the 
effects of gravitational redshifting of an accretion disk 
reflection spectrum, detailed spectral modeling rejects 
this possibility. The addition of a relativistic disk re- 
fl ection component (modeled with the relxill_lp code 
of lGarcfa et al.1120141 ) leads to no significant improvement 
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Figure 5. Data-to-model ratios (left) and (right) resulting from our spectral analysis of the NuSTAR data (FPMA=red; FPMB=blue). 
In all cases, the spectral model includes a two-temperature ICM component with shape (but not normalization) fixed to that found in the 
XMM-Newton analysis. The rows differ in the description of the AGN spectrum. Row-(a) : AGN described as a power-law absorbed by 
a cold column [cABS(PL)]. Row-(b) : AGN described by absorbed power-law plus narrow emission line [cABS(PL+emLINE)]. Row-(c) : 
AGN described by absorbed power-law with a reflection component [cABS(PL+REFL)]. Row-(d) ; AGN described by absorbed power-law 
with reflection, plus the joint absorption/emission signatures of a single-component fast highly ionized wind. Row-(e) : Same as row-(d) 
except for the inclusion of a second ionized absorption component. Error bars are 1-cr. 


in the goodness of fit (Ay 2 = —1 for 5 additional model 
parameters), even if we permit the height of the irradi¬ 
ating source, and the metallicity, ionization state, incli¬ 
nation and inner edge of the accretion disk to all be free 
parameters. 

This iron-K band complexity is not simply due to the 
presence of an iron-edge in reflection. Adding a Comp¬ 
ton reflection component to describe reflection from the 
outer accretion disk or obscuring torus (REFL; see Sec- 
tionUj) improves the fit at 10-30 keV but does not account 
for the spectral complexity in the iron-K band. Indeed, 
since the reflection component forces the iron-Ka line 
to be at 6.4 keV, the overall fit is slightly worse than 
the more phenomenological ICM+cABS[PL+emLINE] 
model. The iron-K band residuals persist even if the ion¬ 
ization state and metallicity of the reflector are allowed 
to be free parameters (model ICM+cABS[PL+REFL] in 
Tabled also see panel-c of Fig. EJ. Since this reflection 
model now includes a Compton hump, the underlying 
power-law continuum is inferred to be steeper (photon 
index T = 1.77 ± 0.03 compared with T = 1.54 ± 0.02 for 
models without the Compton reflection). However, the 
steepening of the primary continuum within this reflec¬ 
tion model now leaves high-energy residuals in the form 
of a hard tail above 40keV (Fig. [SJ panel-c). 

Still, in gross terms, this simple model (ICM emission 


plus an absorbed AGN power-law continuum with asso¬ 
ciated reflection) provides a good description of the NuS¬ 
TAR spectrum. The presence of Compton-hump clearly 
points to the presence of scattering from a Compton- 
thick, or near Compton-thick, structure in the central 
engine of Cygnus A. We return to this issue in Section [6] 

- here we simply point out that this matter must lie 
outside of our line-of-sight since our same spectrum con¬ 
strains the absorption column to the central X-ray source 
to be only Nh = 1.69jj}T x 10 23 cm -2 (approximately 
0.12 Thomson depths). To the best of our knowledge, 
this is the first direct demonstration of Compton-thick 
circumnuclear material in the nucleus of Cygnus A. 

5.2. Modeling the Detailed Iron-K Band Structure : 

Detection of a Fast Ionized Wind 

We now return to the iron-K band complexity noted 
above. Referenced to the fiducial simple reflection model, 
the 4-8 keV spectrum shows correlated sets of residu¬ 
als at the ±5% level, with excursions out to ±10% level 
(Fig. 0 panel-c). The correlated nature of these resid¬ 
uals argues against simple Poisson fluctuations (which, 
for the binning used in Fig. [5j are ±3%). Furthermore, 
these residuals are at a substantially greater level than 
expected effective-area calibration errors which should 
be at the < 1% level based on an analysis of NuSTAR 
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spectra of the C rab neb ul a ma tched to the same off-axis 
angles (also see iMadsen et al.ll2015h . 

Given that Cygnus A has strong ICM emission lines, 
we also need to consider the effects of uncertainties in 
our description of the ICM as well as residuals arising 
from uncertainties in the energy calibration of the FPMs. 
Allowing the shape-parameters of the two-temperature 
ICM model to vary within the 90% error ranges deter¬ 
mined by XMM-Newton has a negligible effect on the fit 
and the residuals (Ay 2 <0.1 for three additional but 
constrained degrees of freedom). If we allow the ICM 
model to be completely free, the normalization of the 
low-temperature component increases by a factor of 100. 
The corresponding impact on the iron-K band, driven 
by the very strong increase of the 6.7 keV line complex, 
has a significant impact on the residuals under discussion 
here. However, this now brings the ICM model into very 
strong conflict with X'MM-Newton, over-predicting the 
soft (< 2keV) emission by a fact of 5-10. We conclude 
that ICM models that are compatible with the XMM- 
Newton data leave these residuals in the NuSTAR fit. To 
address the calibration question, we refit our fiducial re¬ 
flection model permitting a free gain-scale — this can re¬ 
duce the magnitude of the iron-K band residuals if there 
is an energy-scale offset in the iron-K band of 120 eV (in 
both FPM). However, the gain scale of NuSTAR has been 
calibrated on a pixel-by-pixel basis using the onboard ra¬ 
dioactive calibration sources. The expected residual gain 
calib ratio n error this early in the mission is 5E < 20 eV 
(IMadsen et all 12015 ). These considerations suggest that 
the residuals are real and have an astrophysical origin. 

The pattern of residuals suggests that, in addition to 
the spectral features associated with Compton-thick X- 
ray reflection, there are both additional emission and ab¬ 
sorption components. Given that relativistic disk reflec¬ 
tion fails to explain these features, we turn our atten¬ 
tion to wind models. Specifically, we suggest that these 
spectral features correspond to blue-shifted iron-K ab¬ 
sorption together with (redshifted) iron-K emission from 
the back-side of the wind (i.e., an iron-K band P-Cygni 
profile). 

As a first exploration of this wind hypothesis, we add a 
narrow Gaussian emission line (emLINE) and/or a nar¬ 
row absorption line (absLINE) to our base-line model 
ICM+cABS(PL-|-REFL). These three additional model 
fits are presented in Table [2] Adding just the absorption 
line leads to a significant improvement in the goodness 
of fit (Ay 2 = —21 for two additional model parameters) 
but still leaves significant iron-band residuals. The ad¬ 
dition of just the extra emission line produces a much 
greater improvement in the goodness of fit (Ay 2 = — 59 
for two additional model parameters). However, the 
residuals left by this model still indicate the presence 
of absorption, and so it is not surprising the best model 
(ICM+cABS[PL+REFL+emLINE+absLINE]) incorpo¬ 
rates both emission and absorption lines. 

A more physical description of the wind requires 
photoioni zation mod eling. We use the XSTAR code 
([Kallman &; Bautistall2001l) to compute both the X-ray 
absorption and emission spectra from a photo-ionized 
slab of matter (with solar abundances) as a function of 
ionization parameter and column density Nw- Our 
grid of models cover the range £w € (10,10 5 ) erg cm s -1 
in logarithmic steps of A £,w/^w = 0.2 (sufficient to re¬ 


solve the structures in ionization space; [Reynol ds et al.1 
I2012D . and the range Nw g (10 21 ,10 25 ) cm -2 with steps 
A Nw/Nw = 0.2. Guided by the simple-line fitting re¬ 
sults above, we expect velocities in the 0.03 — 0.1c range. 
Thus, we construct the XSTAR models assuming line- 
of-sight velocity spreads of 10,000 kms -1 (formally, this 
is achieved by setting the turbulent velocity parameter 
to 10, 000 kms -1 ). We note that the results described 
below do not change appreciably if we adopt smaller ve¬ 
locity spreads. 

Augmenting our base-line spectral model for Cygnus A 
(ICM+cABSpL+REFL]) with a blueshifted absorption 
component (wABS) and a redshifted emission compo¬ 
nent (wEMIS) of the same ionization parameter £\y pro¬ 
duces a significant improvement in the goodness of fit, 
Ay 2 = —74 for 5 additional model components. The 
iron band residuals are, to a very large degree, described 
by these wind features (Fig. [5] panel-d). Table [2] re¬ 
ports the parameter values for this model. We find 
that, formally, we only obtain a lower-limit on the col¬ 
umn density of the absorber, Nwabs > 3.4 x 10 23 cm -2 , 
and that we require the wind to be highly ionized, 
log£w = 3.2 — 3.9. We note that the XSTAR model, as 
constructed, does not correctly capture Compton scat¬ 
tering and hence does not formally reject Compton-thick 
solutions. Confirming expectations from the simple line 
fitting, the velocity of the absorbing outflow is in the 
range 12, 000 — 29, 000 kms -1 . The emission component 
of the outflow has a fitted column density that is com¬ 
parable with, or maybe slightly less than, the absorption 
component (./Vwemis = 1.0 —10 x 10 23 cm -2 ). This shows 
that the outflow extends a significant solid angle as seen 
from the central X-ray source. 

This description, with a single ionized absorption com¬ 
ponent, still leaves a weak unmodeled absorption feature 
in the residuals (at 7.5 keV observed energy; Fig. [5jl). 
Adding a second ionized absorption component (wABS) 
with its own independent column density, ionization pa¬ 
rameter, and velocity, leads to a further improvement 
in the fit (Ay 2 = —19 for three additional model pa¬ 
rameters). While we have concluded the possible pres¬ 
ence of this second wind component from systematic 
spectral analysis, we note that the presence of (two) 
absorption lines at observed energies of 6.5 keV and 
7.5keV is suggested even by examining the residuals 
from a fit of a simple absorbed power-law to the data 
(Fig. 01). We note however, while it appears to be very 
fast (20,000 — 100, 000 kms -1 ), this second wind com¬ 
ponent has poorly constrained properties. The full spec¬ 
tral model for this two-component ionized wind model is 
shown in Fig.jOj together with a zoom-in on the iron line 
region of the spectrum. 

To a large degree, the poor constraints from the wind 
fit result from strong covariances between model param¬ 
eters. To uncover these covariances, we have performed 
a Goodman-Weare Monte Carlo Markov Chain (MCMC) 
analysis of our two-component wind model as imple¬ 
mented in Jeremy Sanders’ XS PEC-EMCEE pac kage. 
The G oodman-Weare algorithm dForeman-Mackev et al.1 
1201311 simultaneously steps a set of N walkers through 
parameter space, with each walker taking a random step 
in the direction of another (distinct) randomly chosen 
walker. In addition to having excellent convergence prop- 
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Figure 6. Left panel : Our best fitting NuSTAR spectral model for Cygnus A. Shown here are the power-law continuum absorbed 
by both a neutral column and a highly ionized fast outflow (thin black line), the hot ICM component (red line), reflection from ionized 
Compton-thick matter (blue line), and the redshifted re-emission from the ionized wind (green line). Right panel : Folded NuSTAR 
spectrum zoomed in on the iron complex. The FPMA and FPMB data have been co-added and binned for plotting purposes (but not for 
fitting). The solid and dashed lines show the best model fits with and without the fast ionized wind. 


erties, this flavor of MCMC does not require a previously 
defined “proposal” (i.e. probability distributions detail¬ 
ing how each MCMC step should be taken) and is su¬ 
perb for uncovering previously unknown covariances. We 
run 50 walkers for 5000 steps (leading to a 250,000 ele¬ 
ment chain) with a 1000 step burn-in period. Since our 
XSTAR models do not accurately capture the effects of 
Compton scattering, we impose a prior that the column 
density not exceed Nw = 3 x 10 24 cm -2 (i.e. approxi¬ 
mately two Compton depths). While we do not exclude 
the possibility of higher column densities, our current 
models certainly will be a poor description of such winds. 

The results of the MCMC analysis for selected pa¬ 
rameters are shown in Fig. 0 Several interesting co- 
variances are found. We see that the column density 
of the wind components (both absorption and emission) 
are strongly positively correlated with their ionization 
parameter. This is readily understood, at least for the 
absorber, when it is realized that the data are sensitive to 
essentially only one strong absorption line (FeXXVI-Ka). 
For sufficiently high ionization parameters, the fraction 
of iron in that charge state is fFeXXVi cx £,w • an d so 
the column density of that ion is Npexxvi oc Nwabsffw ■ 
Given that observations fix the depth of the absorption 
line and hence the column density of FeXXVI, we expect 
a degeneracy of the form Nwabs oc f,w- Similar argu¬ 
ments drive the degeneracy between Nwemis and £,w- 
Of course, an immediate consequence of these degenera¬ 
cies with is that the two wind column densities are 
also covariant. We note that the wind velocity is free 
from any obvious degeneracies/covariances. 

5.3. Constraints on the High-Energy Continuum Cutoff 

If the X-ray continuum is formed by thermal Comp- 
tonization in an accretion disk corona with electron tem¬ 
perature T e , we expect the power-law spectrum to pos¬ 
sess an exponential cutoff with a characteristic energy 


of E cut « 2 — 3 kT e . Adding an exponential cutoff to the 
continuum component of our best-fitting wind+reflection 
spectral model and refitting produces a slight improve¬ 
ment in the goodness of fit (Ay 2 = —4.1 for one ad¬ 
ditional model parameter). A blind application of the 
F-test would suggest that the cutoff is significant at the 
90% level (but not the 95% level). However, the F-test 
is not well posed in this situation since the null model 
(with no cutoff) is not fully embedded in the expanded 
mode l (with a cutoff), exis ting only in the E cut —> oo 
limit dProtassov et al.ll2002 ). Thus, we conclude that we 
cannot claim the detection of a cutoff. The lower-limit 
on the cutoff energy is E cut > lllkeV at 90% confi¬ 
dence (> 101 keV at 99% confidence). We note that 
Compton-reflection is strongly required even when a con¬ 
tinuum cutoff is included in the spectral model - in¬ 
cluding the cutoff in our single-wind model does not lead 
to a statistically-significant change in the require Comp¬ 
ton reflection fraction (shifting the 90% error range from 
n = 0.25 - 0.51 to n = 0.13 - 0.35). 

5.4. Do We See the Ionized Wind in Previous X-ray 

Data? 

With the 2013 NuSTAR data suggesting the existence 
of a high-column density, high-ionization, and fast wind, 
we must return to the XMM-Neuiton/EPlC-pn dataset 
(§4) and ask whether these wind signatures are also 
present during this 2005 observation. We can also ask 
whether hints of these signatures have been seen in any 
previous X-ray observations. 

Applying the photoionization wind model, we find 
that the 2005 EPIC-pn spectrum requires neither the 
blueshifted absorption nor the redshifted emission from 
a fast wind. With an assumed line-of-sight velocity 
spread of 10, OOOkms -1 and fixing the ionization pa¬ 
rameter and velocity at the best-fitting NuSTAR values, 
the formal limits on the column density of the wind are 
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Figure 7. Results of MCMC analysis of the wind model as applied to the NuSTAR data. The column densities of the absorbing (iVwabs) 
and emitting ( Nwemis ) components of the wind are in units of 10 23 cm -2 . The redshift of the absorbing component (z\y a b s ) is with 
respect to the observer, and so .zwabs = 0.056 correspond to being at rest with respect to the nucleus of Cygnus-A. F is the photon index of 
the primary continuum. The histograms on the left-hand edge of each row show the probability distribution for the parameter in question 
and are normalized to unity. 


Nwabs < 1.2 x 10 23 cm -2 and Nwemis < 4.5xl0 22 cm -2 
(90% confidence), strongly inconsistent with the NuS¬ 
TAR values. 

Thus, we must appeal to time-variability in order to 
reconcile these two datasets. In fact, time-variability 
is clearly required irrespective of the detailed spectral 
modeling — simple PL+emLINE characterization re¬ 
quires the line centroid to shift from 6.39 ± 0.01 keV 
in the XMM-Newton data to 6.34 +qq| keV in NuSTAR 
(a change that exceeds the 20 -30 eV gain calibration of 
NuSTAR ; IMadsen et aLll2015lh 

We suggest that the fast ionized wind has appeared 
during the 7.5 years between these observations. In 
the same period, the 2- lOkeV (de-absorbed) luminos¬ 
ity of the AGN has more than doubled, changing from 
1.7 x 10 44 ergs _1 during the XMM-Newton observation 


to 3.7 x 10 44 ergs 1 during our NuSTAR pointing E 
Over the same time-frame, the cold/neutral intrinsic 
column has dropped from Njj « 3.1 x 10 23 cm -2 to 
Nh ~ 1.4 x 10 23 cm' 2 . Thus, it is tempting to suggest 
that some fraction of the ionized outflowing matter seen 
by NuSTAR originated from the cold material that ob¬ 
scured the AGN during the XMM-Newton observation, 
possibly in (highly non-linear) response to an increase in 
the continuum luminosity. 

What about other X-ray datasets? Suzaku observed 
Cygnus A on 15-Nov-2008 with 46 ks of good on-source 
exposure. The ICM aspects of this dataset have been pre- 

17 The observed/absorbed 2-10 keV AGN luminosity increased 
from 9xl0 43 ergs -1 during the XMM-Newton observation to 2.Ox 
10 44 ergs -1 during the NuSTAR pointing. 
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Observed 

Fluxes 

T 3 —lOkeV = 3.1 x 1CT 11 erg cm -2 s —1 
Tl0-30keV = 5 -6 x to -11 erg cm -2 s —1 
T 30 - 8 O keV = 6.8 X 10 -11 erg cm -2 s —1 

Intrinsic AGN 

U2-i0keV = 3.7 x 10 44 ergs —1 

luminosities 

TlO—30keV = 4.2 X 10 44 ergs -1 


£ 30 —80keV = 5.1 x 10 44 ergs —1 


Table 3 

Observed fluxes (including ICM contributions and effects of 
absorption), and intrinsic AGN luminosities (excluding ICM and 
de-absorbing), based upon our best-fitting model from Tabic [2] 

viously discussed by iSarazin. Finoeuenov &; Wiki d2013t) . 
To examine the core, we have extracted the spectra 
from the three operating X-ray Imaging Spectrometers 
(XISs) on Suzaku using an extraction radius of 3'. Fit¬ 
ting our fiducial ICM+cABS[PO+REFL] model to the 
0.7-10 keV XIS spectra (with ICM parameters allowed 
to be free) we note residuals at the ±10% level in the 6- 
9 keV region. However, there are deviations of the front- 
and back-illuminated XISs at the same level, and there is 
an obvious effect from a background line at 7.5 keV. We 
conclude that the Suzaku-XIS data can neither confirm 
or refute the presence of wind-like spectral residuals. 

The other observatory capable of, in principle, wit¬ 
nessing the spectral signatures of a highly ionized wind 
is Chandra. In the vast majority of Chandra /ACIS data 
on Cygnus A in the archive, studies of the AGN spec¬ 
trum are strongly compromised by severe photon pile-up. 
However, there is a small amount of exposure time (9 ks) 
taken in a short frame-time (0.4 s) mode on 26-May-2000. 
A detailed stu dy o f the AGN using these data has been 
conducted by lYoung et all (120021) . They note the ex¬ 
istence of an anomalous absorption feature blue-wards 
of the 6.4 keV fluorescent emission line which they em¬ 
pirically describe using an absorption edge with a (rest- 
f rame) threshold energy of 7.2 ± 0.1 keV (see Fig. 5 of 
lYoung et all 120021) . While further study of this feature 
was not possible in such a short exposure, it is tempting 
to suggest that this feature resulted from iron-Ka line 
absorption in an ionized wind. We note that, as was 
the case during our NuSTAR observation, the obscuring 
column to the AGN central engine during this Chandra 
observation was lower than that found during the XMM- 
Newton pointing. 

6. DISCUSSION AND CONCLUSIONS 

6.1. Compton Scattering in Circumnuclear Material 

NuSTAR has given us the cleanest view yet of the 
nuclear hard X-ray emission of Cygnus A. The source- 
dominated spectrum across the 3-70 keV band permits 
us to detect spectral curvature consistent with Compton 
scattering from cold, Compton-thick (or near Compton- 
thick) material that may be identifiable with the outer 
regions of the accretion disk or a high column density 
“core” of the obscuring torus of unified AGN schemes. 
In the spectral modeling presented above, the Compton 
scattering and associated iron flu orescence was de scribed 
using the xillver model of ICarcfa et all (j2013P — we 
stress that, while it is usually employed to describe re¬ 
flection from accretion disks, we are using xillver in the 
spirit of a generic, self-consistent, and accurate model of 
X-ray reflection from Compton-thick matter. 

Another view of the Compton scattered compo¬ 


nent can be gained by using the MYTorus model of 
iMurpliv fe Yactoobl (120091) . While possessing less up- 
to-date atomic physics data and a simpler (non self- 
consistent) ionization structure, this model allows us to 
explore scattering/reflecting distributions that are not 
completely Compton-thick. Replacing the xillver com¬ 
ponent with the scattering and iron fine emission de¬ 
scribed by MYTorus in our ultimate, best-fitting two- 
zone ionized wind model, we find that the fit is com¬ 
parably good, with the column density of the scatterer 
constrained to be -/Vn.scat = 1.9+^ x 10 24 cm -2 (0.7- 
2.0 Compton-depths) at 90% confidence, and -Ah, scat = 
1.9l 2 ;f x 10 24 cm -2 (0.5-2.8 Compton-depths) at 99% 
confidence. Taken at face value, this suggests that 
the scattering/reflecting matter is only mildly Compton- 
thick. Confirming our previous expectations, the scatter¬ 
ing/reflecting column density is substantially more than 
that seen in (cold) absorption, demonstrating that the 
scattering matter must be out of the line of sight. 

6.2. The X-ray Loudness of Cygnus A 

We use our best fitting (wind±reflection) model to de¬ 
rive intrinsic AGN luminosities, correcting for the con¬ 
tribution from the ICM and de-absorbing. The results 
are shown in Tabled] The 2-10keV X-ray luminosity is 
3.7 x 10 44 ergs -1 , al most 10% of the bolo metric luminos¬ 
ity as calculated by iPrivon et al. (120121 ) on the basis of 
the radio-to-IR SED. Unl ess t he I R-calo rimetry assump¬ 
tion underlying the IPrivon et al.1 (12012D work fails due 
to some unusual geometry, we conclude that Cygnus A 
is X-ray loud, with a 2-10 keV bolometric correction 
of only 10. However, given the fact that Cygnus A 
likely has an Eddington ratio of only 10 -2 , this behav¬ 
ior falls in linejwith t hat f ound in other broad-line AGN 
(iVasudevan &; Fabianl 120091) . 

6.3. Physical Properties of the Wind 

Our most novel finding in the NuSTAR spectrum is 
the presence of subtle features indicative of a high- 
column density (Nwabs > 3 x 10 23 cm -2 ), highly ionized, 
(£ ~ 2500ergcms -1 ), fast (v = —2.0tQg x 10 4 kms -1 ) 
wind. The fact that we see possible red-shifted emission 
from this wind suggests that it is a wide-angle outflow, 
subtending a significant solid-angle of the sky as seen 
from the central X-ray source. Here, we discuss the in¬ 
ferred physical properties of this wind. 

Suppose that the wind signatures originate from clouds 
with characteristic hydrogen number density n and vol¬ 
ume filling factor / at a characteristic radius ro- The 
observed column density is N\y = nrof and the ioniza¬ 
tion parameter is £ = Li/nr g = L,/ /N^rg (Z^ is the 
ionizing luminosity of the central source). The mass flux 
is the wind is then given by 

M\y = Llfm p pLrQVwn ( 1 ) 

= Llm P iJ.NwroVw , (2) 

where Q is the solid angle subtended by the wind as 
viewed from the center of the system, vw is the velocity 
of the wind, and p is the average particle mass of the 
wind per proton (in units of the proton mass m p ). We 
now make the assumption that the wind is flowing at its 
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local escape velocity so that the characteristic radius is 



With this, we can write useful expressions for the mass 
flux 

M w = 2 Q,m p <?nr g N w Vw , (4) 

scalar momentum flux (Pw), 

Pw = M\yvw = 2flm p c 2 /j,r g Nw, (5) 

and kinetic energy flux ( Lk ), 

Lk = ^M w vw = flrripC 2 nr g N w vw, ( 6 ) 

of the wind. 

Putting our lower limit on column density Nwabs > 
3.4 x 10 23 cm -2 into eqn. Ogives 

P w > 1.4 x 10 36 dyne « 10 , (7) 

where we have used a black hole mass of M = 2.5 x 
10 9 Mq and assumed that the wind covers a solid an¬ 
gle of fl = 7T. Thus, the momentum flux in the wind 
exceeds the photon momentum flux by an order of mag¬ 
nitude. Including the information on the wind velocity 
also gives us lower bounds on the mass and energy flux 
(using eqns. Q] and [6]), 

M w > 7.7 M 0 yr - 1 » 110 , ( 8 ) 

L k > 1.7 x 10 45 ergs -1 « 0.42L bo i- (9) 

We can immediately see that the mass flux in the wind 
exceeds the accretion rate onto the black hole unless the 
radiative efficiency is 77 < 9 x 10 -3 , and that the energy 
flux is at least 42% of the bolometric luminosity. These 
bounds are based on the assumption that the observed 
wind velocity is the local escape velocity (thereby allow¬ 
ing us to localize the wind in radius). If, instead, we 
made the weaker assumption that the wind is moving 
at no less than its local escape velocity, the momentum, 
mass and energy fluxes further increase. 

Given these lower limits on wind fluxes, we can 
make some comment about acceleration mechanisms. 
For radiative-driving, the fact that Pw > 10Pboi/c 
would require acceleration of a very Compton-thick wind, 
7 7 , > 10 , th at co mpletely surrounds the radiation source 
(|Revnoldsll2012f) . While not ruled out by these data, this 
would require a revision of our basic picture for the cen¬ 
tral engine structure of AGN. Magneto-centrifugal accel¬ 
eration is possible provided that the inward mass flux in 
the accretion disk exceeds the mass flux in the wind (in 
order for the disk to generate sufficient torque to generate 
the wind). To make this compatible with our measured 
limits on the mass flux requires a low accretion efficiency, 
77 < 0.01, as might be expected for these low Eddington 
ratios due to the formation of an optically-thin advective 
accretion flow. If Cygnus A does indeed possess an ad¬ 
vective accretion flow, we also expect some contribution 
to the wind acceleration to c ome from thermal driving 
ijBlandford fc Begelmanlll999t) . 


Of course, Cygnus A also possesses very powerful 
relativ istic jets. Usin g the dynamics of the radio co¬ 
coon, llto et al.l (1200811 estimate a total jet power in 
the range Lj « (0.7 — 4) x 10 46 ergs _1 ; this exceeds 
the radiative luminosity by a factor of 2-10. What¬ 
ever process dri ves these relativistic jets, whether it be 
accreti on dBlandford k, Pavnel 119821 1 or black hole spin 
dBlandford fc ZnaiekHl977U we need only channel < 25% 
of the power into a wide-angle sub-relativistic outflow in 
order to explain the observed wind. 

6.4. Implications for AGN Feedback 

A large body of prior wo rk has i denti fied two dis¬ 
tinct modes of AGN feedback (lFabianll2012H . Relativistic 
radio-jets from AGN in massive ellipticals and the bright¬ 
est cluster galaxies in cooling core clusters appear to pro¬ 
vide “maintenance mode” feedback, preven tin g a coo ling 
catas trophe in their hot atmospheres dPeterson fe Fabianl 
1200611 . On the other hand, powerful sub-relativistic 
winds from luminous quasars can sweep molecular gas 
out of a gala xy and quench star formation following a 
major merger faopkins et al.ll2006tlTombesi et al.ll2015f ). 
These are normally considered to be mutually exclusive 
forms of feedback. This notion is partially motivated by 
drawing an analogy between AGN and stellar mass black 
hole X-ray binaries (BHBs) in which the mutually exclu¬ 
sive occurrence o f disk winds and relativistic jets is q uite 
well established (|Ponti et al.ll2012t King et al1l2013f) . 

Cygnus A appears to break this mold. We clearly see 
powerful jets interacting with the ICM of the cooling 
core cluster. However, we now also see a powerful and 
wide-angle sub-relativistic wind with sufficient energy 
and momentum to sign ificantly impact the host galaxy 
(Hopk ins fe Elvisll2010t) . We should not be too surprised 
that some breakdown in the AGN-BHB analogy exists. 
Important aspects of the accretion physics do not remain 
invariant as we scale the mass up from a stellar-mass 
to a supermassive object. AGN accretion disks are sig¬ 
nificantly more radiation-pressure dominated than their 
BHB cousins, raising the possibility that their stability 
properties are different. Additionally, being cooler, AGN 
disks have much higher opacities at their pho tosp heres 
raisin g the possibility of line-driven winds (lLaor &; Davisl 
120141) . Line-driving may prime the fast wind that we see, 
lofting material off the surface of the disk where it can be 
strongly photo-ionized by the central X-ray source and 
magneto-centrifugally accelerated. 
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